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We investigate the microscopic origins of nonlinear radio frequency (RF) response in supercon- 
ducting electromagnetic resonators. Strong non-linearity appearing in the transmission spectra at 
high input powers manifests itself through the emergence of jump-like features near the resonant 
frequency which evolve towards lower quality- factor with higher insertion loss as RF input power is 
increased. We directly relate these characteristics to the dynamics of localized normal regions (hot 
spots) caused by microscopic features in the superconducting material making up the resonator. A 
clear observation of hot spot formation inside a Nb thin film self-resonant structure is presented 
by employing the microwave laser scanning microscope, and a direct link between microscopic and 
macroscopic manifestations of nonlinearity is established. 



PACS numbers: 74.81.-g, 74.25.N-, 74.62.Dh 
I. INTRODUCTION 



Microwave resonators utilizing superconducting thin 
films have enabled development of various nonlinear de- 
vices including bifurcation amplifiers 1 , hot electron mix- 
ers^, and kinetic inductance photon detectors^ among 
other examples. They have been crucial elements in 
circuits designed for quantum computation^. Recently, 
these on-chip resonators have started to appear in meta- 
material designs promising new possibilities such as nega- 
tive refractive index^, electromagnetically induced trans- 
parency 7 , tunable THz response^, and parametric ampli- 
fication of negative index photons 9 . 

Nonlinearity is manifested in many forms, including 
a reduction in quality factor, an increase in insertion 
loss, shift of the resonant frequency^, or appearance of 
abrupt jumps in the transmission/reflection character- 
istics^. Understanding the microscopic origin of such a 
response is essential. Various mechanisms have been pro- 
posed to explain the nonlinear characteristics at elevated 
microwave power, including both intrinsic and extrinsic 
types. Among them, hot spot formation is found to be 
one of the major mechanisms, which was both theoreti- 
cally and experimentally studied in superconducting strip 
conductors and resonators^. 

A clear and compelling connection between the mi- 
croscopic models and the macroscopic picture of nonlin- 
earity revealed by microwave measurements is still lack- 
ing. Here we establish this connection through a study of 
the radio frequency (RF) response of a two-dimensional 
spiral resonator made up of 200 nm Nb thin film on 
a quartz substrate via the Laser Scanning Microscopy 
(LSM) technique. 



II. EXPERIMENT 

The resonant structures have been designed as engi- 
neered meta- atoms (planar spirals) of RF metamateri- 
als intended to create a strong magnetic response in 
the sub- 100 MHz rang o 13 i 14 . Prior results show that 
the spirals act as very compact self-resonant structures, 
supporting up to 10 half- wavelength standing waves of 
RF current^. The free space stimulation of these spi- 
ral resonators is done by sandwiching them between two 
RF antennas made up of loop-terminated coaxial cables 
whose other ends are connected to an RF network ana- 
lyzer. Exciting the spirals results in Lorentzian type res- 
onant peaks in transmission for low levels of RF stimulus 
power—. The peaks do not change in shape or position 
over a broad range of input powers. When the spirals are 
strongly driven, the resonant characteristics begin to de- 
velop notches near the center of the transmission peaksii. 
Those nonlinear features evolve with RF power and are 
believed to be associated with the existence of non super- 
conducting regions in the Nb thin films causing resistive 
hot spot formation. 

The high-power properties of superconducting strip 
resonators are very sensitive to the geometrical structure 
of the current-carrying conductors^. Our continuous cir- 
cular spirals have a unique geometry (see the left inset of 
Fig. 1) in which the currents flowing in neighboring turns 
are parallel, and approximately equal in magnitude, at 
least for the first few resonant modes, in contrast with 
the current flow profile in meander-line resonators used 
for photon counting 4 . In the spiral, the magnetic fields 
generated by flowing currents largely cancel in the region 
between the strips, leading to a self- field pattern nearly 
parallel to the plane of the spiral. This geometrical ef- 
fect greatly reduces the accumulation of RF currents at 
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the edges of the strip. Moreover, the continuous arc of 
the spiral allows for a relatively homogenous and uni- 
form flow along the windings compared to sharp-cornered 
structures that have a tendency to pile up the currents 
at the corners and edge o 16 i 17 . 



III. RESULTS AND DISCUSSION 
A. Global Microwave Transmission Measurements 

The main panel of Fig. 1 shows the global microwave 
transmission response of a Nb spiral resonator for a set 
of RF input powers at 4.4 K. The spiral has an outer 
diameter, D Q of 6 mm and 40 turns; the strip width, w 
and spacing, s between the turns are 10 fim. The funda- 
mental resonance gives a pronounced transmission peak 
at ~ 75.84 MHz for low and moderate RF powers (see 
Fig. 2 in Ref^). Up to -2 dBm, the resonant peaks al- 
most overlap, i.e., transmission, |S2i(/)| is a Lorentzian 
function of frequency and does not show any significant 
difference in shape with increasing power (black curves 
in the right inset). Those curves falling into the linear 
response regime correspond to a state where Nb is in the 
hotspot-free superconducting phase. At a critical input 
power (-1.5 dBm in this case), |6 , 2i(/)| makes a sharp 
transition from one Lorentzian curve onto another with 
higher insertion loss and lower quality factor, Q, as fre- 
quency is scanned near resonance (the highest red curve 
in the right inset). 

As the driving frequency is swept from lower values 
towards the resonance frequency, fo, the circulating RF 
power builds up in the resonator and the circulating cur- 
rents I C irc m the spiral will increase in magnitude, as 
I max l^l + [{f-h)/5ff where Imax is the max- 
imum current induced in the resonator and Sf is the 
3 dB bandwidth. For low/medium RF input powers, 
no instantaneous power value causes excess dissipation 
because even the maximum current at resonance is low 
enough to maintain the linear characteristics. For the 
high input power case, the cavity is over-driven in the 
vicinity of fo and thermal effects cause the breakdown of 
superconductivity. The over-stimulation of the resonator 
launches a new dissipation mechanism and the resonant 
response moves onto a different curve with lower Q corre- 
sponding to initial stages of a hot spot. With increasing 
input power this transition occurs at progressively lower 
frequencies where the dissipative mechanism is activated 
(the red curves in the right inset of Fig. 1). At still higher 
powers, a second discrete dissipative mechanism is ac- 
tivated, leading to another set of discontinuous curves 
(blue curves). At even more elevated powers one has 
either a large number of discrete dissipation centers ac- 
tivated, or an expansion of the original hot spots, giving 
rise to a flattening of the transmission spectrum (orange 
curves), and eventually to a dip or crater (brown curves). 

The fo of the Nb spiral remains roughly constant with 
increasing power, unlike what is seen in temperature de- 
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FIG. 1: (Color online) Global microwave transmission re- 
sponse from a D — 6 mm spiral with 40 turns and w =s — 10 
\im for a set of RF input powers at T = 4.4 K. Right inset: A 
closer look at the linear and nonlinear regimes evolving with 
increasing input power. Black Lorentzian type curves corre- 
sponding to power values from -21 dBm to -2 dBm overlap, 
progressively evolving into the red curves for power values be- 
tween -1.5 dBm and +0.6 dBm then blue curves for powers 
from +0.7 to +1.1 dBm, orange curves from +1.2 dBm to +2 
dBm and brown curves from +2.2 dBm to +3.6 dBm. Top 
Left inset: An optical image of a spiral resonator with D = 
6 mm and 40 turns. Bottom Left inset: Zoom in the area 
marked in the top left inset showing the Nb spiral windings 
in detail. Dark strips are Nb turns in the picture. 



pendence data where the resonance shifts to smaller fre- 
quencies 7 . At the extreme case of input power, the dis- 
sipation is so large that superconductivity is quenched 
and the entire resonant feature is wiped away from the 
spectrum (see the data at +27 dBm). Note that upon re- 
ducing the power the resonant features are re-established. 



B. Laser Scanning Microscopy 

To develop a microscopic understanding of the nonlin- 
ear resonant features in the transmission data, we exam- 
ined the spiral resonators with the LSM techniqu e - 17 ! 18 . 
For these measurements, the Nb spiral is stimulated by 
RF signals well below its transition temperature, T c and 
is scanned by a focused laser probe. The spiral resonator 
is supported on a sapphire disk along with a thermometer 
and heater located on a metallic platform nearby, inside 
an optical cryostat. The RF excitation has been done 
by two RF loops in a manner identical to that described 
above. A diode laser beam of 5 mW power and a wave- 
length of 645 nm illuminates the spiral surface. For the 
large-scale LSM images (up to 10x10 mm 2 scanned area 
with 500x500 points), a collimated laser beam is focused 
onto a spot with a diameter of 20 \im that produces an 
estimated temperature difference, 5T of 1.2 K at most 
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relative to an area on the sample unexposed to the laser 
beam. For more detailed examination of the spiral sur- 
face (within a maximum scanned area of 250 x 250 /im 2 ), 
a 1.6 jim diameter laser beam is used. The raster scan- 
ning of the spot is done with 0.2 jim steps, much smaller 
than the laser spot diameter. However, the amplitude of 
the laser is modulated at a frequency of 100 kHz, which is 
significantly smaller than the inverse phonon escape time 
(on the order of nanoseconds), thus creating a steady 
state temperature profile^. Moreover, the beam inten- 
sity is reduced by neutral-density optical filters down to 
1 mW/mm 2 , so the temperature rise is reduced to ST= 
0.1 K. 

The LSM technique basically images photoresponse 
(PR) which comes from the changes induced by pe- 
riodically modulated laser power in the transmission 
characteristics (~(d\S2i(f)\/dT)5T) of the spiral, syn- 
chronously detected with a lock-in amplifier. The re- 
sponse can be decomposed into inductive and resistive 
components^. Since the temperature is well below the 
T c of Nb for the measurements shown here, the mod- 
ulation of penetration depth A (thus the inductive re- 
sponse) is small, so the contrast in the LSM images 
is mainly due to resistive PR, which is proportional to 
S(f R s \ 2 Jft F (x,y)dS), where the integral is carried out 
over the area of laser-induced perturbation. The PR is 
therefore a convolution of changes in resistance below T c , 
5R S weighted by the local value of RF current density 
squared, J^ iF (x,y). 

Large-scale (7x7 mm 2 ), low temperature (4.8 K) LSM 
images of the Nb spiral presented in Fig. 1 are shown in 
Fig. 2(a) and (b). The excitation of the spiral is done at 
the fundamental frequency of ~ 73.6 MHz^ at which the 
first emerging hot spot is detected at +10 dBm (a). By 
the time the power is ramped up to +20 dBm, the entire 
sample is heated and the resistive response allows imag- 
ing of a convolution of SR S and the current distribution in 
the entire spiral (b), showing a single half- wavelength of 
standing wave current spanning the length of the spiral. 
Most of the PR comes from the middle windings, and di- 
minishes near the inner and outer edgesii The brighter 
areas observed on the left side and the lower half of the 
spiral in (b) are caused by the structural/morphological 
defects with enhanced SR S . 

The upper half of the spiral is free of physical imper- 
fections (no defect is visible with 0.5 fim resolution re- 
flected light microscope), thus we have concentrated on 
a small portion of that area to examine the appearance 
of the first microscopic defect. Figs. 2(c)-(f) are a de- 
tailed examination of the 30x30 fim 2 area on the Nb spi- 
ral marked in (a) and show the evolution of the first hot 
spot with increasing RF power at 4.8 K. At +10 dBm, 
we observe two bright spots in their initial phase at the 
opposite edges of a single strip in the spiral (c), rem- 
iniscent of vortex-antivortex formation at the edges of 
the current carrying superconducting stripes 2 ^. As the 
power is increased, the isolated spots expand and start 
to merge together (d) where presumably phase-slip lines 
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FIG. 2: (Color online) LSM PR images of the nucleation 
and RF power evolution of a hot spot, (a) Detection of the 
first microscopic defect causing the hot spot formation on 
the Nb spiral resonator whose details are discussed in Fig. 1 
at +10 dBm, 4.8 K. (b) The hot spot shown in (a) heats 
the neighboring area and creates a resistive response coming 
from the entire spiral at +20 dBm. The length scale bar in 
(a) is the same for (b). (c)-(f) Zoom in of the 30x30 /am 2 
area marked in (a) at progressively higher input powers. The 
color-scale bar in (d) applies to all LSM images shown in Fig. 
2, the maximum LSM PR corresponds to 0.3346 mV for (a), 
0.1857 mV for (b), 3.6 \iV for (c), 11.4 yV for (d), 23.2 fiV 
for (e) and 12.4 jiV for (f). The length scale bar in (c) applies 
to the images in (c)-(f). 



(i.e. isothermal nonequilibrium regions)^! are formed. 
Further increase in RF power causes a change of the ar- 
rangement of these phase-slip lines (e) due to the redistri- 
bution of current flow along the strip. This follows from 
the nucleation of a normal domain across the strip inter- 
facing with resistive superconducting domains [the bright 
regions established along the wings of a butterfly pattern 
seen in (e)] 22 . With even further increase in power the 
normal domain expands along the strip and finally, the 
entire strip shows a resistive response (f). 

The formation of a phase-slip line in the DC resistive 
state is signaled by the emergence of step-like features in 
the /(V) characteristics. Although the phase-slip lines 
discussed in this paper are generated by the combination 
of RF excitation and laser power, we see a very similar 
picture. The data presented in Fig. 1 with gradually de- 
grading Q-factor with increasing RF input power make 
a good analogy to the characteristics of the DC bias case 



with changing differential resistance as the voltage is in- 
creased. 

We further studied the structure of the hot spot with 
another sub- 100 MHz Nb spiral at 4.5 K. The spiral ex- 
amined in Fig. 3 has 45 turns, D Q = 5 mm and w =s =10 
fim. Here, the laser beam passes through the hot spot on 
a 20 jim line scan [shown with 'X' in Fig. 2(c)]. The fre- 
quency is swept through resonance (/o~ 88-2 MHz) along 
the X-line cut, to create a position vs. frequency image 
of the hot spot at a fixed RF input power. Figs. 3(a)- 
(d) show RF power evolution of the hot spot within a 
20 fim line as frequency is scanned in a narrow band 
centered on /q, while the corresponding PR profiles at 
a fixed position, xq near the center of the hot spot are 
shown in (e)-(h) for both directions of frequency sweep. 
At dBm, in the initial stages of the hot spot, one sees a 
strong Lorentzian PR{f ', xo) along with a linear |5 f 2i(/)| 2 
[(a) and (e)], reminiscent of that at low power shown in 
Fig. 1. In the case of a spatially invariant 5R Sl the dis- 
sipative term of PR(f) originates from the temperature 
dependence of J^ F (f)R s (T)^ meaning that it can be de- 
scribed by a simple difference between two |6 f 2i(/)| char- 
acteristics 5\S 2 i(f)\~Jl F (f)[tts(T 2 ) -i? a (Ti)], where T x 
is the temperature of the non-radiated area while T 2 is 
the temperature of the area radiated by the laser probe 
on the Nb film. At low RF powers, the temperature ef- 
fect of the laser on transmission will be small, hence the 
PR traces out a Lorentzian-like curve, as seen in (e). 

As the power increases, the hot-spot produces a non- 
linear behavior of 5R s (Jrf,T) causing peculiar jumps 
in \S 2 i(f)\ whose trend is mimicked by PR(f). At +4 
dBm [(b) and (f)], in the middle of the hot spot, one 
first sees a Lorentzian-like response at the onset of the 
frequency sweep, but then a breakdown occurs at a crit- 
ical frequency, and a crater develops in the PR(f) trace 
similar to that seen in the higher power \S 2 \(f)\ char- 
acteristics in Fig. 1. Note that these line cuts are not 
hysteretic; i. e. forward (blue curves) or backward (red 
curves) sweep of the frequency reproduce the same char- 
acteristics. 

An important question arises: why does the LSM PR 
drop in the crater? Is it due to a decrease in SR S or Jrf, 
or both? In fact, both mechanisms should contribute 
since LSM PR ~(dR s /dT)(5T)+(dJ RF /dT)(5T). The 
dissipative component of LSM PR is proportional 
to both SR S and changes in the normal (non- 
superconducting) component of Jrf, SJrf- The latter is 
reduced in the resistive state due to the diminished super- 
fluid density. In this case, a crater-like pattern is formed 
due to time averaging of LSM PR coming from the local- 
ized regions where resistive superconducting and normal 
domains are side by side, as in phase-slip processes 21 i 23 . 
When a fully normal domain is formed, SR S ~ 5Rn, 
where SRn is the very small change in normal state re- 
sistance with temperature, thus the LSM PR drops sig- 
nificantly, as seen in (g). 

A higher input power (+5 dBm) LSM imaging of the 
same x-line scan is shown in (c). Note that there are 
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FIG. 3: (Color online) Frequency and space dependence of 
LSM PR at (a) dBm, (b) +4 dBm, (c) +5 dBm (d) +6 
dBm on a Nb spiral with D a = 5 mm, 45 turns, w =s = 10 
jim over a 20 jim x-scan. The frequency is swept between 86.5 
and 89.5 MHz through the fundamental resonant frequency, 
fo = 88.2 MHz. The images in (a)-(c) show only forward fre- 
quency sweeps whose direction is indicated with blue arrows. 
In (d), different shade color-maps point out both direction 
sweeps, shown with blue (forward) and red (backward) ar- 
rows. Corresponding PR(f,x) profiles at a fixed position, 
x — xo of the laser probe are shown in (e)-(h). The blue 
dots are the data collected through forward frequency sweep 
whereas the red dots are those through backward sweep. The 
|*S f 2i(/)| 2 at low power, dBm is shown in (e). The color- 
scale bar in Fig. 2(d) applies to all LSM images shown in 
Fig. 3 as well, the maximum LSM PR corresponds to (for- 
ward/back sweep)-(1.61/1.62) mV for (a); (1.49/1.5) mV for 
(b); (1.54/1.55) mV for (c); and (1.12/1.4) mV for (d). 



now two shoulders and a central valley at resonance in 
the PR profile shown in (g). Our interpretation is that 
there is a RF critical (still partially superconducting) 
state with significant 5R S for excitation frequencies on 
the shoulders. In the PR valley, the material has become 
normal, mimicking what we have observed in Fig. 2(e). 

With more power, at +6 dBm, the central valley at 
resonance becomes deeper (d) and at some point the PR 
drops to zero (h). Though away from the resonance some 
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FIG. 4: (Color online) Evolution of global transmission with 
small increments in RF input power (a), Corresponding LSM 
PR averaged over 100 fim line scans across the center of a 
hotspot. Note that |5 f 2i(/)| is measured simultaneously with 
LSM PR(f) when frequency is changed very slowly. A nearly 
exact correlation between |&i(/)| and PR data is evident. 



recovery is observed, heating is so significant that back 
and forth sweeps of frequency can not reproduce the same 
PR profile. Note the hysteresis in (h); the direction of 
the frequency scan is indicated by blue and red arrows. 
The size of the resistive domain expands along the x-axis 
gradually with increasing power and reaches its maxi- 
mum at +6 dBm (d). 

Line cuts at different x locations in Fig. 3 show a con- 
tinuous variation of LSM PR(f) from Lorentzian-like, to 
flat-top, to cratered, which follows what is seen in the 
global |5 f 2i(/)| measurements shown in Fig. 1. These 
data provide very clear observation of such a unique con- 
nection between local temperature change and global 
RF-power dependent properties of a superconducting 
resonant system. 

To further support the claims made above, Fig. 4 shows 
the evolution of global |6 , 2i(/)| characteristics (a) and 



averaged 100 jim line scans of LSM PR (b)-(d) with in- 
creasing power at 4.5 K, on a 40-turn, D Q = 6 mm spiral. 
The linear |S , 2i(/)| characteristics occurring at +10 dBm 
[black curve in (a)] turns into nonlinear curves with small 
increments in RF power. The curves are shifted with 
equal offset in the vertical direction for clarity. With 
increasing input power, the LSM PR appears over pro- 
gressively broader ranges of frequency around resonance 
[(b)-(d)]. The development of PR contrast at the hot 
spot corresponds almost exactly to the deviations of the 
global |6 , 2i(/)| response from the linear/Lorentzian be- 
havior. This demonstrates quite clearly the direct cor- 
respondence between hot spot formation and nonlinear 
response in superconducting resonators. 

IV. CONCLUSIONS 

In conclusion, we presented experimental evidence 
linking the global nonlinear microwave transmission char- 
acteristics of a superconducting resonator with hot spots 
induced inside the material making up the resonator. 
The hot spots are imaged with LSM technique at various 
stages of development and their influence on the global 
transmission characteristics is documented. 
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